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ABSTRACT 

We report the discovery of two binary millisecond pulsars in the core-collapsed globular cluster M30 using 
the Green Bank Telescope (GBT) at 20 cm. PSR J2 140-23 lOA (M30A) is an eclipsing 1 1-ms pulsar in a 4-hr 
circular orbit and PSR J2140-23B (M30B) is a 13-ms pulsar in an as yet undetermined but most likely highly 
eccentric (e > 0.5) and relativistic orbit. Timing observations of M30A with a 20-month baseline have provided 
precise determinations of the pulsar's position (within 4" of the optical centroid of the cluster), and spin and 
orbital parameters, which constrain the mass of the companion star to be m2 > 0.1 M©. The position of M30A 
is coincident with a possible thermal X-ray point source found in archival Chandra data which is most likely 
due to emission from hot polar caps on the neutron star In addition, there is a faint (V555 ~ 23.8) star visible in 
archival HST F555 W data that may be the companion to the pulsar Eclipses of the pulsed radio emission from 
M30A by the ionized wind from the compact companion star show a frequency dependent duration (oc z^"" 
with a ^ 0.4-0.5) and delay the pulse arrival times near eclipse ingress and egress by up to 2-3 ms. Future 
observations of M30 may allow both the measurement of post-Keplerian orbital parameters from M30B and 
the detection of new pulsars due to the effects of strong diffractive scintillation. 

Subject headings: galaxy: globular clusters: individual: M30 — pulsars: individual: PSR J2140-2310A — 
radio continuum: stars 



L INTRODUCTION 

Globular clusters (GCs) produce millisecond pulsars 
(MSPs) at a rate per unit mass that is up to an order-of- 
magnitude greater than the Galaxy(e.g. Kulkarni & Anderso3 
Qj96). Due to the relatively large distances of GCs (several to 
tens of kiloparsecs), the low intrinsic luminosities of MSPs, 
and the fact that most MSPs are members of compact binary 
systems, the discovery of new cluster pulsars requires long 
observations with the largest radio telescopes and computa- 
tionally intensive data analysis. 

The discovery of new cluster pulsars is interesting be- 
cause of the wide variety of science that can result from 
using them as sensitive probes into the natures of the pul- 
sars themselves and the clusters in which they live. Re- 
cently, cluster pulsars have been used to probe properties of 
GCs, such as th e mass-to-light ratios in cluster cores (e.g. 
iFreire et'al]l200 1a: D'Amico et al. 2002), cluster proper mo- 
tion (|Freire et al. 2003), and the ionized gas content in 47 Tu- 
canaelFreire et al. 2001b). For binary pulsars, timing obser- 
vations have measured relativistic effects such as the advance 
of peria stron (and th erefore the total mass) in the 47 Tuc H 
system jFreire et"ani2b0^ . and probed the companion winds 
and eclipse mechanisms for several known eclipsing MSPs 
(e.g.lp'Amicoetal. 2001b; Possenti et al. 2003). The pre- 
cise astrometry provided by MSP tuning has allowed the op- 
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tical ide ntification of severaljjinary MSP companions (e.g. 
iFerraro et al. 2001 ; Edmonds et al. 2002), which is crucial for 
determining the nature of the companion sta rs, and the X- 
ray id entification of many MSP systems (e.g. iGrindlav et all 
l2002h . which gives us useful information on pulsar emission 
and neutron star cooling mechanisms. Finally, many theo- 
rists have predicted that truly exotic pulsar sys tems, such as a 
pulsar-black hole binary (e.g.lSi gurdssonl2003i) . will be found 
in GCs. 

After a flurry of GC pulsar discoveries in the 1980s 
and early 1990s, the number of known cluster pulsars re- 
mai ned virtually constant (~35) until 2000 (for a review, 
see iKulkarni & AndersonI Il996h . Over the past several 
years, however, the art of searching for radio pulsars in 
GCs has undergone a renaissance due to the development 
of very sensitive (i.e. low-noise and high-bandwidth) 20- 
cm receivers (e.g. Stavelev-Smith et al. 1996) and the in- 
creasing availability of the high performance computing re- 
sources required to conduct sensitive but specialized searches 
for binary millisecond pulsars in observations with du- 
rations of several hours (e.g. Eo hnston & Kulkarni I991I 
iRansom. Eikenberrv. & Middleditch 2002). The Parkes radio 
telescope has been particularly productive as of late with the 
discovery of at least 24 millis econd pulsars in 8 GCs, most 
of which are in binaries ( Camilo et al. '"2OOO'; D'Amico et alJ 
2001a b; Possenti et al. 2001; Ransom 2001; D'Amico et al] 
2002; Lorimer etalJ2 003: Possenti et al. 2003). 

In the past two years, the recently upgraded Arecibo 
telescope and the new 100-m Green Bank Telescope 
(GBT) have become available, and several new algo- 
rithms have been developed to improve s earch sensitivi- 
ties to bina ry MSPs in compact orbits (Chandler 20031 
iRansom. Cord es. & Eikenberrv 2003). Using the GBT and 
one of these new techniques, Jacobv et al. (2002) have re- 
cently reported the discovery of three new binary MSPs in 
M62 (see also Chandler 2003). With these advances in mind, 
we undertook a major survey of rich and/or nearby GCs at 



2 



Ransom et al. 



20 cm that are visible from Arecibo and/or the GBT. We are 
analyzing the data using modern search algorithms and tech- 
niques. This is the first in a series of papers describing the re- 
sults from these observations. Specifically, we focus on "First 
Science" data from the GBT, taken in the Fall of 2001 toward 
12 GCs. Additional discoveries — including at least one other 
new MSP in Ml 3 discovered in the GBT data described below 
— and results from the rest of the project will be presented 
elsewhere**. 

2. SEARCH OBSERVATIONS AND DATA ANALYSIS 

In 2001 September and October, we used the Gregorian fo- 
cus 20-cm receiver (1.15-1.73GHz usable bandwidth with 
Tsy, ^ 25 K) on the GBT to observe 12 different clusters for 
either 4hrs (M2, M4, M75, M80, M92, and NGC 6342) or 
8hrs (M3, M13, M15, M30, M79, and Pall) each. Since the 
GBT beam has a FWHM of nearly 9' at 20 cm, single point- 
ings fully covered each cluster Samples from two orthogo- 
nal polarizations were transmitted through fiber optic cables 
to the observatory control room where they were fed to ei- 
ther one or two Be rkeley-Caltech Pulsar Machines^ (BCPMs; 
lBackeretalJl997h . The BCPMs are analog/digital filterbanks 
which 4-bit sample each of 2x96 channels at flexible sam- 
pling rates and channel bandwidths and which can sum the 
two polarizations in hardware if requested. 

We observed each cluster using 96x1. 4 MHz channels of 
2 summed polarizations centered at 1375 MHz and sampled 
every 50 fis. For those observations using a second BCPM 
(which included the observations of M30), we recorded 
96x1.0 MHz channels with summed polarizations sampled 
every 36 /is at a center frequency of 1490MHz. In total, these 
initial search observations generated ~0.5 TB of data, which 
were stored on DLT IV magnetic tapes. 

Reduction of the data took place on a dedicated Beowulf- 
style cluster of 52 dual-processor 1.4 GHz AMD Athlon 
workstations located at McGill University. The available 
computing power allowed us to attempt extensive interference 
removal (see below) and to conduct fully coherent accelera- 
tion searches'" for binary pulsars in all clusters. In previous 
GC pulsar projects acceleration searches were only used on 
clusters where pulsars — and therefore the dispersion mea- 
sure (DM) to the cluster — were known a priori, due to the 
computational expense of searching over the additional "ac- 
celeration" parameter. 

Since these data were taken as part of the GBT "First Sci- 
ence" program, the telescope was still in the early stages of 
commissioning and the radio frequency interference (RFI) en- 
vironment was quite bad. Each of the observations contained 
large quantities of both persistent and transient broadband and 
narrowband interference. The very strong Lynchburg, Vir- 
ginia airport radar, with a rotational period of approximately 
12 s, was particularly destructive and effectively eliminated 
any sensitivity we should have had to weak pulsars with pe- 
riods longer than ^ 100 ms. Due to the challenges involved 
in excising this RFI, approximately one quarter of these data 
remain to be fully analyzed. 

We searched the raw data from each cluster for RFI in both 
the time and frequency domains as a function of both BCPM 

^ An up -to-date catalog of GC pulsars can be found at 
[htt^^^^ww . naic . edu/ $ \ sim$pf reire/GCpsr . html 
^"^l^^^^^^www . gb . nrao . edu / $ \ sim$db acker 

Acceleration searches correct either a time series or its Fourier transform 
to account for large apparent per iod derivatives in the si gnals from binary 
pulsars in compact orbits (see e.g. ljohnston & Kulkamll9$B . 



channel and time (in units of ^^10 s). Sections of data contain- 
ing fluctuations significant at the 4-a level in the frequency 
domain or lO-cr level in the time domain were masked using 
the running median of the appropriate BCPM channel. In ad- 
dition, we computed and Fourier transformed non-dispersed 
topocentric time series in order to identify strong and obvious 
terrestrial interference (i.e. occurring at "un-natural" frequen- 
cies such as 11.1 Hz) which we ignored in subsequent stages 
of the analysis. 

After applying the interference masks, we de-dispersed the 
data over a range of DMs from ^-^50-200% of the predicted 
D M of the cluster (usi ng the galactic electron density model 
of lTavlor & Cordesll993iD . or, for clusters with known pulsars, 
from ~90-110% of the average DM of those pulsars. The 
clusters that we observed with the GBT were chosen partly 
because their known or predicted DMs were J^lOOpccm""*, 
thereby keeping the dispersive smearing across each BCPM 
channel to <0.5 ms at 20 cm. The stepsize in DM was chosen 
to produce no more than ~0. 2-0. 3 ms of dispersive smear- 
ing across the 134.4 MHz bandwidth. The final stage of data 
preparation involved removing samples from or adding sam- 
ples to the time series as appropriate in order to account 
for the observatory's motion with respect to the solar system 
barycenter (using the DE200 ephemeris of Standish 1982). 

After Fourier transforming each de-dispersed and barycen- 
tered time series and removing the previously identified RFI 
signals and their harmonics, we searched the data for pul- 
sations in three steps: 1) using acceleration searches of the 
full-length observation to maximize our sensitivity to isolated 
or long-period binary pulsars, 2) with acceleration searches 
of overlapping '--^20 min and ^60 min segments of the obser- 
vation to improve our sensitivity to pulsars in compact bina- 
ries or which exhibit significant diffractive scintillation, and 
3) using phase-modulation (or sideband) searc hes for pulsars 
in ultra-compact (Port ^S2hr) binary systems jRansom et alJ 
2003). For the acceleration searches, we used a modern 
Fourier-domain matched-filtering code that includes Fourier 
interpolation and the incoherent harmonic summi ng of the 
first 1 , 2, 4, and 8 harmonics of any potential signal jRansorrJ 
The acceleration search code allowed the highest 
summed harmonic to linearly drift by up to 170 Fourier 
frequency bins during the portion of the observation being 
searched (lower harmonics drift a smaller fraction of bins 
equal to their harmonic number over that of the the highest 
harmonic). We folded the raw data for interesting accelera- 
tion search candidates over a range of nearby periods, period 
derivatives, and DMs in order to maximize the signal-to-noise 
ratio. For phase-modulation candidates, we used Fourier- 
domain matched-filte ring of orbital templates as described in 
iRansom et al .1(120031) . 

Given the highly variable RFI environment, the often ex- 
treme effects of diffractive scintillation (see ^4.2> . and the 
sensitivity losses caused by uncorrected orbital motion for bi- 
nary pulsars, it is difficult to place a hard limit on our over- 
all search sensitivity. However, based on search simulations, 
our detections of known cluster pulsars, and the radiometer 
equation, we estimate that we were sensitive to "normal" mil- 
lisecond pulsars with flux densities in the range ^-^50-100 /iJy. 
Our sensitivities to sub-millisecond pulsars and ultra-compact 
binary MSPs are factors of ^2-3 worse than the above values 
due primarily to the effects of dispersive smearing across the 
relatively wide BCPM channels for the former and the inco- 
herent nature of the phase-modulation search technique for 
the latter. For slow pulsars (P > 0.1 s), our sensitivities were 
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Fig. 1. — Pulse profiles from the 7.8-hr discovery observation of the 11- 
ms binary pulsar PSR J2 140-231 OA (M30A, left) and the 13-ms eccentric 
binary pulsar M30B (right). The Doppler effects of each orbit have been 
removed, and two complete profiles are plotted for clarity. The data were 
taken with the GBT on 2001 September 9 using the BCPMl machine (see 
O at a center frequency of 1370 MHz with an effective time resolution of 
~0. 12 ms. The greyscale below the average pulse profiles shows the consis- 
tency of the pulsed emission as well as some effects due to interference as 
a function of time. Strong scintillation and portions of at least two eclipses 
of M30A are evident, including dispersive delays to the pulse arrival times at 
ecHpse ingress and egress (see i|4.3l . 

degraded even more due to the devastating effects of RFI, and 
in particular, the Lynchburg radar 

3. DISCOVERIES AND TIMING OBSERVATIONS 

One of the most promising clusters in our search was 
M30 (NGC 7099; l = 2T.\i,b = -46°. 84), due its "post- 
core-collapsed" nature and relatively high central luminos- 
ity density. M30 has a measured core of <2.5" as deter- 
mined by HST ( Yannv et al. 1994), and an age and distance 
of 12.3 Gyr and D=9.0±0.5kpc (appare nt distance modulus 
m-My = 14.90 ± 0.05, £(B-y) = 0.039: ICarretta et all2000h 
respectively as determined from Hipparcos data. Significant 
evidence has been found of mass segregation within the clus- 
ter ( Guhathakurta et al. 1998). 

While performing acceleration searches of ^20min seg- 
ments of BCPMl data for M30 taken 2001 September 9, we 
discovered two highly variable and highly accelerating pul- 
sar candidates near 1 1 ms (M30A; see |3 and 13 ms (M30B; 
see whose signal-to-noise ratios peaked at DMs very near 
25 pc cm"^. Upon examining the BCPM2 data from the same 
time period, we confirmed both of these candidates as pulsars 
with very high significance. Soon thereafter, we determined 
"local" timing solutions for these data consisting of a compact 
4.2-hr orbit for M30A and a high-order polynomial expansion 
in spin frequency for M30B which showed each pulsar visible 
for ^90% of the 7.8 hr observation (see Fig.^. 

Our discovery of two MSPs within the 9' beam of the GBT 
at 20 cm — at least one of which is within 4" of the cluster 
center (uMio = 21^ 40™ 22?16 ± 0"2, Smm) = -23° 10' 47."6 ± 
0."2; J2000; see j4!4l — and which have measured DMs 
('-^25pccm~^' ) similar to that predicted for the cl uster us- 
ing both the iTavlor & CordesI (^23 dc cm"^: 11993ft and the 



NE2001 (^41 DC cm"^: ICordes & Lazioll2002ft models of the 
galactic electron density, firmly establish these as the first 
known pulsar members of M30. 

In 2002 March, we began a series of monthly 4-8 hr ob- 
servations of M30 using the GBT/BCPMl setup described in 
^ We have been unable to use the second BCPM machine 
during these sessions due to problems with the instrument. 
The observations were intended to establish phase coherent 
timing solutions for both pulsars and to allow us to search 
for new pulsars in M30 whose measured flux densities were 
temporarily boosted above our detection limits via scintilla- 
tion. On several occasions, we observed M30 with 48 MHz of 
bandwidth centered at either 575 MHz or 820 MHz in order to 
determine the DM of the pulsars more accurately, take advan- 
tage of a better RFI environment (particularly for 820 MHz), 
and to search for pulsars with spectral indices too steep to be 
detected near 1400 MHz. 

4. PULSAR J2140-2310A (M30A) 

For pulsar M30A we folded the raw data from the timing 
observations modulo the predicted pulse period given our best 
orbital ephemeris and accumulated the resulting pulse pro- 
file until enough signal was available to allow an accurate 
measurement of the phase of the pulsation. We measured 
the pulse ph ase by cross-correlating (in the Fourier domain; 
lTavloill992:') the observed profiles with a high signal-to-noise 
template profile determined during a period of scintillation- 
heightened flux density. Absolute times-of-arrival (TOAs) 
resulted by referencing the phase measurement to the start 
time of each observation as recorded from the observatory 
clock. We later used TEMPO" to correct these times to 
the UTC(NIST) time standard with data from the Global Po- 
sitioning System and to transform them to the solar-system 
barycenter using the DE200 planetary ephemeris (StandisQ 
Il982l) . 

Using TEMPO, we then fit the TOAs iteratively to a model 
incorporating the pulsar position, DM, spin period (P), period 
derivative (P), and the Keplerian orbital parameters for a cir- 
cular orbit; the projected semi-major axis (x = asin//c), the 
orbital period (Port), and the time of the ascending node (Tasc)- 
Attempts to fit for eccentricity resulted in a 95% confidence 
upper limit of e < 1.2 x 10""^. The timing model produced 
RMS residuals of ^24 /is (see Fig.|2} and resulted in the pa- 
rameters shown in Tabled 

4.1. Cluster Accelerations 

It has been well established (e.g. lPhinnevlll992h that at 
least four different effects can contribute at the few percent 
or greater level to the measured spin-down rate, P, of a MSP 
in a GC. Typically parameterized in terms of "accelerations," 
and ignoring the possible effects of a nearby star or planet 
perturbing the system, the measured acceleration of a pulsar 
is 

P~"P c c ~' 

where Po is the intrinsic spin-down of the pulsar, fl£,Gc/c 
and a£,Gai/c are the line-of-sight accelerations caused by the 
gravitational potentials of the cluster and the Galaxy, and 
«pm/c = n^D/c is the apparent acceleration caused by the 
transverse Doppler effect where fi is the pulsar's measured 

'^ |http: / /pulsar . Princeton. edu/tempol 
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Fig. 2. — (Top) Pulse phase residuals for pulsar M30A as determined after 
fitting for the timing model shown in Table [Tl using TEMPO. All measured 
TOAs occurring between orbital phases 0.12-0.38 were excluded from the 
fit in order to minimize systematic effects caused by the pulse delays during 
ecHpse ingress/egress. The RMS residual for the 408 TOAs is 23.7 fis. (Bot- 
tom) All measured phase residuals plotted as a function of orbital phase. For 
a circular orbit, ecHpses are expected to occur during supeiior conjunction at 
a phase of 0.25. Eclipse dela ys of up to several ms are evident during eclipse 
ingress and egress (see i|4.3l . 



Table 1. Parameters for PSR J2140-2310A 



Parameter 


Value 


Right Ascension, a (J2000) 


2ih 4ora 22?40610(46) 


Declination, 5 (J2000) 


-23° 10' 48"7936(97) 


Dispersion Measure (pc cm"'' ) 


25.0640(41) 


Pulsar Period, P (ms) 


11.0193290688805(67) 


Period Derivative, P (s/s) 


-5.181(20)x 10-20 


Epoch (MJD) 


52162.0 


Orbital Period, Pg,-!, (days) 


0.17398746418(34) 


Projected Semi-Major Axis, x (It-s) 


0.2349416(48) 


Eccentricity, e 


<1.2xl0-^ 


Epoch of Ascending Node, Tasc (MJD) 


52161.94552243(68) 


Span of Timing Data (MJD) 


52162-52850 


Number of TOAs 


408 


Weighted RMS Timing Residual (ps) 


23.7 


Flux Density at 1400 MHz, 5i4oo (mJy) 


0.08(3) 


Derived Parameters 


Mass Function, fi (Mq) 


0.000459967(28) 


Minimum Companion Mass, m2 (Mq) 


>0.10 


Radio Luminosity, L1400 (mJykpc^) 


6.5(2.5) 



Note. — Numbers in parentheses represent the formal 2a uncertainties 
in the last digit as determined by TEMPO, after scaling the TOA uncertainties 
such that xt/'^-^- The values for fi and m2 were derived assuming a pulsar 
mass of 1.4 M0. 



proper motion and D is its distance JShklovskiill 1 9701) . Typ- 
ically, the intrinsic spin-down and cluster acceleration terms 
are of the same order, while the Galactic and proper motion 
terms contribute at the 10% level or less. 

In the case of M30 and M30A, we know that P/P = 
-4.7 X 10"'^ s"' and can estimate a£,Gai/c ~-8.5 x 10"'^ s"' 
assuming a spherically symmetric Galaxy with a flat ro- 
tation curve (Phinney 1993) and upm/c — 1-3 x 10""^ s"' 
given the 7.8 mas yr"' of proper motion measured for M30 



bv lDinescu. Girard. & van Altenal jl999ll) . To an accuracy of 
- 10%. lPhinnevUl992lll993l) showed that 



max - 



|fl,,Gc| l.lGE«e±) 



(2) 



where !](< 8j^) is the projected surface mass density within 
the radial position of the pulsar in the cluster, 3.6" for 

M30A). If the projected surface luminosity density is avail- 
able from optical measurements, then we can use Equation|2l 
to place a lower limit to the projected mass-to-light ratio of 
the cluster within the radius of the pulsar position. 

Using the power-law relation for the total-li ght V-band sur- 
face flux density in the core of M30 from Guhathakurt a et all 
{1998), we estimate the surface luminosity density within the 
radius of the position of M30A to be ~ 2.0 x IO'^LqPC"^. 
Substituting and rearranging Equation ^ implies that M/L > 
0.51 Mf7i/Lf7i in th e core of M30, which is typical for GCs. 

iPhinnevI j 19921 Il993h also showed that to within 10% if 



flf.GC 



3v^(0) 



2cD 



(3) 



where Qc is the core radius of the cluster (which we have as- 
sume d to be ^^1.8" based on the range quoted bv lYannv et alJ 
fT99l . and v^O) is the one-dimensional velocity disper- 
sion for the cluster core (9.4±2.5kms"' for M30 from 
iGebhardt et al. 1995). Since the measured P/P is negative, 
the intrinsic "acceleration" must be 



Po \ae,Gc\ 
— < max — ■ — 
P c 



P 
+ -- 
P 



«PM_ (4) 

c 



If we account for the l-cr error bars in D, Vi(0), and the 10% 
accuracy of Equation |3] for M30A, Po/P < 1.4 x 10""' s"'. 
Similarly, the surface magnetic field strength of M30A is 
Bs = 3.2 X 10'.''(PPo)'/2G < 4.2 x 10'*G, the characteristic 
age Tc =P/(2Po) > 1.1 X 10* yr, and the spin-down luminos- 
ity E = 47r^IP„/P^ < 4.6 X lO^'* erg s ' assuming the canonical 
/ = 10"^^ g cm^. Each of these values is typical for MSPs. 

4.2. Pulsed Flux 

The monitoring observat ions have rev ealed significant 
diffractive scintillation (e.g. lRickettlll97'i) which results in 
large-scale variations of the measured flux density from 
M30A over time scales of one to several hours, and over 
bandwidths of 50-100 MHz. By integrating the pulsed signal 
above the average off-pulse levels and comparing the mea- 
sured noise variance with that described by the radiometer 
equation, we have estimated the 20-cm flux density (S'woo) 
to an accuracy of ^30% during 73 1-hour sub-integrations 
which occurred during un-eclipsed portions of the pulsar's 
orbit. We did not detect the pulsar during four of the sub- 
integrations, so for those epochs we assumed a flux den- 
sity of 1/2 that of the weakest definite detection we have 
for the pulsar ('-- ^ 16 /xJy), similar to the procedure used by 
Camil o et alJ (I2OOO). The measurements, displayed as a cu- 
mulative plot in Fig. 13 show the expected exponential distri- 
bution for strong diffractive scintillation (e.g. Rickett 1977ft 
with an average flux density of ~0.08(3)mJy. We made sim- 
ilar measurements for 8 1-hour intervals at 575 MHz and 10 
1-hour intervals at 820 MHz and found flux densities of ap- 
proximately 0.13(4)mJy and 0.12(4)mJy respectively. These 
measurements correspond to a rather flat radio spectral index 
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Fig. 3. — A cumulative plot of 73 measurements of the 20-cm flux den- 
sity of pulsar M30A. Each measurement comprised an hour of integration 
time with the BCPM/GBT setup as described in ^ The thick curve shows 
a cumulative exponential distribution with an average of 87 fi]y. The mini- 
mum detectable flux density in 1 hr was ~16/xJy. During six of these 1-hr 
integrations that occurred outside of predicted eclipse times, the pulsar was 
undetected. For these points, we assumed a flux density of half the minimum 
detectable value, or 8 fily. 



for this pulsar of -0.6;!^ y. With such a low average flux den- 
sity, it is not surprising that M30A (and the less-luminous 
M30B; s ee © was not detected in earlier imaging surveys 
jHamilton. Helfand. & Becker 1985, using the VLA) or pul- 
sation searches (Biggs &Lvn&,199 6, using the 76-m Lovell 
Telescope at Jodrell Bank) of the cluster 

At the distance of M30, the 20-cm flux density for M30A 
corresponds to a radio luminosity L1400 = 6.5(2.5) mJykpc-. 
Comparing this number to 40 MSPs in the most recent ATNF 
pulsar catalog (Manchester et al., in prep)'^ with measure- 
ments of L1400, 11 catalog MSPs with measurements of 
L400 and assuming a typical spectral index for each of -1.6 
(^orimeret al. 1995), and the 14 MSPs in 47 Tucanae with 
estimates of L1400 in Camilo et al. (2000), we find that M30A 
is in the most luminous ^^30% of known MSPs, and more than 
twice as luminous as the median L1400 of ~2.7 mJy kpc^. 

The lower frequency observations of M30A have revealed 
that the amplitude of the interpulse (as measured with respect 
to the main pulse) increases with decreasing radio frequency. 
While the separation in phase at each of the three observing 
bands 1400/820/575 MHz seems to be constant at --42%, the 
ratio of the flux contained in the main pulse to that in the in- 
terpulse is approximately 1.9/1.2/1.0 at 1400/820/575 MHz, 
respectively. Such frequency dependent evolution of pulse 
component amplitudes — while the pulse phases of the com- 
ponents remain constant — seems to be a com mon occurrence 
among pulsars (see e.g. iKramer et alJll999l and references 
therein). 

4.3. Eclipse Properties 

M30A exhibits consistent and total eclipses at 1400, 820, 
and 575 MHz for approximately 20% of its orbital period 
centered on superior conjunction (i.e. orbital phase 0.25, see 
Figs.^andlJJ- The presence of ecUpses impUes that the orbit 

'^ |http : / /www, atnf ■ csiro ■ au/ research/pulsar /psr cat 7| 



is significantly inclined (/ > 30°) and therefore the mass of the 
companion star is almost certainly in the range 0.10 < mi < 
0.21 Mq. As such, M30A (f = 1 1 .0 ms, P^rb = 0. 174d) is very 
similar to the three other eclipsing MSPs with 0. 1-0.25 Mq 
companions and orbital periods of a few hours, all of which 
are located in GCs: PSR B1744 -24A in Terzan 5 (P = 
1 1 .6 ms, P„rb = 0.076 d: lLvne et alJl9 90). PSR J 0024-7204W 
in 47 Tucanae (P = 2.35 ms, P^rb = 0.133 d; ICamilo et alJ 
l2000l) . and PSR J1701-30 06B in M62 (P = 3.59ms, P^rb = 
0.145d; Possenti et al. 2003). 

Scintillation and pulse delays at eclipse ingress and egress 
(see below) make the determination of the exact duration of 
eclipses as a function of observing frequency difficult. From 
a single good 820 MHz observation, we measure the eclipse 
duration to be 0.23(2) Po,;,, while for four good observations 
near 1400MHz, the eclipse duration is 0.l8{l)Porb- This 
corresponds to a frequency dependence of the eclipse dura- 
tion with a ^ 0.4-0.5, very similar to that measured by 
iFruchter et all (11990) for PSR B 1957-1-20. Such a value is also 
consistent with the cyclotron absorption based eclipse model 
of KhechinashviU, MeUkidze, & Gil (2000), where the ab- 
sorption occurs in the magnetosphere of the degenerate com- 
panion star which gets continuously injected with relativistic 
particles from the pulsar wind. 

During times of increased signal-to-noise ratio due to scin- 
tillation, we observe pulse phase delays of up to 2-3 ms 
during eclipse ingress and egress which imply the pres- 
ence of an additional electron column density of A^^. 7.4 x 
lO'^A, tnst'Qfj^ cm"^ in the eclipse region, where A, tns is the 
pulse delay in ms, and z/ghz is the observing frequency in 
GHz (see Fig. |2ji. For our observations of M30A, this cor- 
responds to A'f > 2 X lO'^cm"^. Since the orbital separation 
of the system is a ~ I.4R0 (for / > 60°) an eclipse dura- 
tion of ^20% of the orbital period corresponds to a physi- 
cal size of the eclipsing region Re of '^O.9R0. The addi- 
tional electron density near the eclipse boundaries is there- 
fore > 3 X 10^ cm"-' which corresponds to a plasma density 
of p£ > 5 X 10"'^ g cm"-* if fully ionized. 

As in the cases of the other eclipsing MSPs with low- 
mass companions, the physical sizes of plausible companion 
types is much smaller than Re- Assuming / = 60° and there- 
fore m2.60° = O.I2M0, hydrogen and helium white dwarfs 
(WDs) have radii of -008 Rg or -^0.03 Rq respectively 
( Shaoiro & Teukolsk \IT98l . while a zero-age main sequence 
star would have a radius of ~O.13R0. Even the Roche lobe 
radius Re — 0.26 R© (e.g. Eggleton 1983) is several times 
smaller than Re, implying that the eclipses must be due to 
wind- or magnetos nhere-related a ctivity at several times the 
companion radius (IFruchter et alj [r990: Stappers et al. 199^ 
iKhechinashvili et alJl200(1l) ." If this wind is emitted isotrop- 
ically, the plasma density at the eclipse boundaries corre- 
sponds to a mass loss rate M = AttR^ePev^, where is the 
velocity of the ionized wind from the companion. If we as- 
sume vv — 10^cms"\ which is the order of the escape ve- 
locity from the surface of a presumed WD companion, the 
corresponding mass loss rate is M 4 x 10"'^ Mq yr"'. As 
has been found with the other eclipsing systems, unless the 
ionized fraction of the eclipsing wind is small (< 0.4), the 
companion star will not be ablated in a Hubble time. 

The nature of the eclipse mechanism for M30A is currently 
difficult to constrain given that we know neither the true spin- 
down luminosity of the pulsar nor the nature of the companion 
star. Additionally, the low average flux density of M30A has 
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allowed only a few good measurements of the eclipse ingress 
and egress. However, it does seem unlikely that the eclipses 
are caused by dispersive smearing of the pulses, since the ex- 
cess Ne that we measure near the eclipse boundaries would 
cause <1 ms of smearing across the bandwidth of our obser- 
vations and would therefore not smear the pulse enough for 
it to appear eclipsed. While our initial measurement of the 
frequency dependence of the eclipse duration seems to be in 
reasonable agreemen t with the cyclotron absorption model of 
iKhechinas hvili et better constraints on the eclipse 

mechanism will have to await higher sensitivity observations 
as well as good measurements of the eclipse properties at sev- 
eral other observing frequencies. 

4.4. Optical Observations 

We have attempted to identify the companion of M30A in 
archival HSTPNFPC2 observations of M30. We have used the 
observations in the F336W, F439W and F555W filters (here- 
after f/336, B439 and Vsss) of GO-5324 (1994 March 1) and the 
observations in the C/336, V555 and F814W (hereafter /814) fil- 
ters of GO-7379 (1999 May 31 and June 1). The exposure 
times of these observations are 200 s in f/336 , 80 s in B439 and 
16 s in V555 for the GO-5324 program and 11 600 s in f/336, 
1 192 s in Vg55 and 1676 s in /814 for GO-7379. These images 
were reduced and photometered with the HSTphot 1 . 1 pack- 
age bv Dolpbia C20 00), following the recommended proce- 
dures in the HSTphot manual'-'. 

The HST/WFPCl images were placed onto the Interna- 
tional Celestial Reference System (ICRS) for direct compari- 
son against the position of the pulsar (Tabled)- A detailed de- 
scription of this procedure is presented in Bassa et al. (2,00^. 
In short, a 8' x 8' section of a 4 minute V-band image, taken 
2000 August 28 with the Wide Field Imager ( WFI) at the ESO 
2.2 m telescope at La Silla, was placed onto the ICRS using 
58 stars from the USNO CCD Asti'ograph Catalog (UCAC; 
IZacharias et a l. 2000.). We fitted for zero-point position, scale 
and position angle, and the astrometric solution has rms resid- 
uals of 0.06" in both right ascension and declination. This 
solution was transferred to the WFPC2 images by matching 
stars on the WFI with those on the WFPC2 chips, where we 
corrected the WFPC2 pixel positions for geometric distortion 
an d placed them on a maste r frame following the prescription 
of lAnderson &King lll2()0l. A total of 99 stars were used for 
the GO-7379 dataset (105 for GO-5324), giving rms resid- 
uals of 0.05" to 0.06" in both coordinates. The final (la) 
uncertainty in the astrometric solution of the WFPC2 images 
is about 0.07" to 0.08" in right ascension and declination. A 
4" X 4" section of an image comprising 1044 s of the V555 
data from GO-7379 is shown in Figure |4] The image shows 
a single star within the 0.25" (95% confidence) error circle 
centered on the M30A position. The star is offset from the 
pulsar position by -0.09" in right ascension and -0.07" in 
declination and is a plausible companion to M30A. 

HSTphot identified and measured the magnitude of the 
candidate companion to M30A in 13 of the 48 individual V555 
images from GO-7379, but failed to detect it in the U^if, and 
/814 images. The S/N ratio of the V555 detections is very low, 
roughly between 2-5, and hence the uncertainties in the mag- 
nitude determinations are large (typically 0.4 magnitudes). 
The average V555 magnitude of the candidate from the 13 de- 
tections is 23.8 ±0.1. We folded the magnitude measure- 
ments at the orbital period of the pulsar using the ephemeris 

'^ |http : / /www.noao ■ edu/ staff /dolphin /hstphot/| 




Fig. 4. — Optical lAs^-band (lett) and X-ray (right) images of portions of 
M30. North is towards the top in each image and east to the left. (Left) A 
4" X 4" HST/WFPC2 F555W (V555) image with a total exposure of 1044 s 
centered on the radio position of M30A. The error circle has a radius of 0.25" 
and is the 95% confidence region for the position of M30A in the optical 
frame. A single star is present in the error circle, right and below of the 
center The countour levels correspond to 50, 100, 200, 500, 1000, and 2000 
counts. See ^3 for more details. (Right) A 50ks, 32" X 32" C/icmdra 
ACIS-S3 0.2-3 keV image of the central portion of M30. Light grey pixels 
indicate a single detected photon, while darker pixels indicate two or more 
counts. Contour levels correspond to 5, 15, 50, and 150 counts. The error 
circle has a radius of 0.6" which is the approximate error in the absolute 
astrometry of the Chandra image. The astrometrically corrected position of 
the optical center of the cluster as determined by G uhathakur ta et al. 1 199^ 
is marked with a black cross. We consider it likely that the ^5 eve nts j ust 
to the northwest of the nominal pulsar position are from M30A. See i|4.5l for 
more details. 



in Tabled] but found no significant variability, possibly due 
to the large uncertainties in the magnitudes. We note that the 
non-detections occur roughly uniformly in orbital phase and 
therefore cannot be explained exclusively by orbital phase de- 
pendent variability. These measurements are consistent with 
visual inspection of the candidate position on the individual 
V555 images. The non-detections of the candidate companion 
in f/336 and /gw place the object on or bluewards of the cluster 
main sequence. 

Whether or not this candidate is the actual companion to 
M30A, it is unlikely that the true compani on is similar t o 
that of PSR J1740-5340A in NGC6397 ( Ferraro et alJ200ir) . 
That star is as luminous as stars at the main sequence turn 
off, while the U^^d magnitude of the candidiate companion 
to M30A is at least 5 magnitudes fainter than that. If the 
companion is a main sequence star, it is mor e likely com- 
pai-able to PSR J0024-7204W in 47 Tucanae JCamilo et alJ 
2000), whose binary companion was recently identified in 
HST observations on the basis of large-amplitude variabil- 
ity which rnatched the orbital period and phase of the pulsar 
(Edm onds et alJ li2002). Since M30A has similar orbital char- 
acteristics to 47 Tuc W, it is plausible that the companion to 
M30A shows optical variability as well. The current V555 data, 
though, rule out variability at the 1 .6 mag level as seen from 
47 Tuc W. A future series of HST/ ACS images might unam- 
biguously identify our candidate as the companion to M30A 
on the basis of lower amplitude variability or detect a more 
likely counterpart. 

We have followed the recommendation of 
Guhathakurta et al. (1998) and matched their astrometry 
against our astrometric solution of the GO-5324 dataset. We 
find that their reference star (No. 361 1) has an absolute posi- 
tion of a2000 = 21^ 40" 22?314 a nd Sjmn = -23° 10^ 40"10 . 
We have furthermore matched the Guhathakurta et all J1998I) 
positions of the 40 stars in their Table 1 against our positions 
and solved for zero-point offset, scale and rotation. The 
[Guhathakurta et al.. (1998) frame had to be shifted by -1 .25" 
and 2.27" in R.A. and in Deck, scaled by a factor of 0.988 
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in both axes and rotated by -0.76° around the new position 
of the reference star, to result in our astrometric frame. The 
rms residuals of the transformation were about 0.006" in 
both coordinates. In o ur absolute astrometric frame the 
iGuhathakurta et al.l (|1998) position of the M30 centi'oid has 
^2000 = 21*^ 40" 22?16 and ^2000 = -23° 10' 47f'6. We use this 
position for the M30 cluster center throughout the paper. 

4.5. X-ray Observations 

We have also attempted to identify pulsar M30A in a 50 ks 
Chandra ACIS-S3 observation centered on the core of M30 
and taken 2001 November 19 (OBSID 2679). We used the 
CIAO''' software package (v3.01 with CALDB v2.23) to ap- 
ply the most up-to-date aspect, charge transfer inefficiency, 
and gain map corrections to produce a new Level=2 events 
file from the archival Level=l events. From that file we kept 
all events with energies in the range 0.3-3 keV in order to 
produce the 32" x 32" image of the cluster core shown in Fig- 
ure |4] Just to the north-west of the nominal pulsar position 
there appears to be a weak (^5 counts with energy <3keV) 
source. While this "source" was not detected by the WAVDE- 
TECT algorithm, given th at several eclipsing MSPs have been 
identified in X-rays (e.g. Grindlav et al. 2002: 'Stappe rs et aP 
12003), and that the timing position of M30A is within the 
~0."6 absolute astrometric error circle for Chandra, it is plau- 
sible that the observed X-ray flux originates from M30A. 

If we assume a blackbody m odel (using PIMMS' ^) with 
a single kT 0.22 keV (as in Grindlav et al. 2002) and a 
Nh for M30 of 3.5 x lO^^'cm"- (Dickev & Lockman 1990), 
the five measured events in the 0.3-3 keV band correspond 
to an equivalent unabsorbed flux from the source of ~ 5 x 
10"'^ ergs cm"^ s"', which at the distance of M30 corresponds 
to a luminosity of ~ 5 x lO'^'ergss"'. This is comparable to 
the X-ray luminosities found for the MSPs in 47 Tucanae and 
NGC 6397 which are believed to be due to emission from hot 
polar caps on the NSs (G rindlav et al.lfe002 ). Alternatively, 
the emission could be due to thermal bremsstrahlung (TB) 
in the hot ionized wind of the companion star where the ra- 
dio eclipses take place. While a TB model with T I keV 
yields an almost identical luminosity of 6 x 10-"' ergs s~', 
the implied emission measure would be EM '--^ 2 x lO^'' cm""*. 
If we assume the plasma to be fully ionized and uniformly dis- 
tributed in a spherical region of diameter Re, the total volume 
of the plasma is 1 .3 x 10^^ cm^ and the corresponding elec- 
tron density is 1 x 10'' cm"-'. This value is ^ 3000 times 
larger than the measured electron density near the eclipse 
boundaries ( ^4.3> . Since the eclipsing material is being con- 
stantly replenished (based on the fact that size of the eclipsing 
region is much larger than even the Roche lobe radius of the 
companion) and since we measure no dispersive delays out- 
side of the orbital phases corresponding to the eclipse region, 
the very large electron densities implied by TB effectively rule 
it out as the source of the X-ray emission. An additional ar- 
gument against the TB model comes from the fact that the 
isolated MSPs in 47 Tucanae show X-ray emission yet have 
no companion stars to generate the required hot ionized wind. 

5. PULSAR J2140-23B (M30B) 

For each of the timing observations, we searched the data 
using the same techniques described in ^both for new pul- 
sars and for a re-appearance of M30B. In addition, we per- 

http : // asc ■ harvard, edu/ciaol 
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Time (MJD - 52162) 

Fig. 5. — The predicted time dependence of tlie measured spin period of 
M30B given four representative Keplerian orbital solutions from the family 
shown in Fig.|^ The thick black line shows the measured spin period during 
the discovery observations of 2001 September 9. 



formed limited folding searches of 20-min to 1-hr data seg- 
ments for M30B over a range of periods centered on its dis- 
covery spin period and a range of plausible period derivatives 
that would be caused by the accelerations of the possible or- 
bits discussed in ^5.11 In no case have we either definitively 
re-detected M30B or confirmed any additional pulsars in the 
cluster 

Our lack of subsequent detections of pulsar M30B can be 
explained in one of two ways. The most likely explanation is 
simply that M30B has a significantly lower average flux den- 
sity than M30A, and that its discovery was the result of an 
extreme diffractive scintillation event. Given the exponential 
distribution of diffractive scintillation events, however, and 
the implied low probability of observing extreme scintillation 
events, it is unlikely that M30B has a flux density more than a 
factor of r^2 below that of M30A. This implies that if diffrac- 
tive scintillation is the cause for our non-detections in blind 
searches, then the pulsar would be visible in much of our data 
if we could search for it optimally (i.e. folding on the correct 
orbital ephemeris). We note that M30A is not detectable via 
blind searches in approximately half of our observations. 

The other possibility for the non-detections of M30B is 
unusual eclipses. It could be that M30B has an intrinsic 
luminosity similar to or even greater than that of M30A, 
but its companion causes irregular and possibly long-term 
"eclipses" of the pulsed signal in the same manner as the 
GC MSPs PSR B1744-24A in Terzan 5 (Lvneetal. 199^ 
Nice & Thorsett 19921 and PSR J1740-5340 in NGC 6397 
(.D'Amico et al.-200l5) . 

5.1. Orbit Constraints 

Without post-discovery detections of M30B we are cur- 
rently unable to unambiguously determine the orbital pa- 
rameters of the system. However, examination of the spin- 
frequency behavior during the more than seven hours that the 
pulsar was visible in the 2001 September 9 discovery obser- 
vation immediately led us to conclude that we could not fit a 
sinusoid to the frequency behavior (see Fig.|5}, and hence that 
the orbit was significantly eccentric. 

We measured 26 TOAs from the BCPMl and BCPM2 dis- 
covery observations and used TEMPO to fit a timing model 
consisting of the DM of M30B and a Taylor expansion of 
the measured pulsation frequency in time. Reasonable solu- 
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tions providing "white" timing residuals required fits to the 
initial spin frequency (/) and at least the first four frequency 
derivatives (/, /, f^^\ and The best 4-derivative solu- 

tion, at an MJD(UTC) epoch of 52161.997395, and a DM of 
25.09±0.12pc cm-\ is 

/ = 77.003746(14)Hz, 

/ = -8.557(59) X IQ-^Hzs"', 

/ = 5.68(16)x 10-"Hzs-^ 
/3) =-4.06(25) X 10"'^Hzs-^ 
/^' = 1.79(19) X 10"'''Hzs"^, 

and provided RMS residuals of ~53 /is and the pulse profile 
shown in Fig.^ 

Given the timing solution above, w e applied the or bit- 
inversion technique described by Joshi & Rasid ([l9^) to 
convert the polynomial-based solution into a family of Ke- 
plerian orbital solutions as a function of orbital eccentricity, e 
(assuming that the intrinsic / for the pulsar is negligible dur- 
ing the observation). A lack of physical solutions from the 
technique provided a rough lower limit to M30B's eccentric- 
ity of e > 0.45. Using the inversion-based family of solutions 
as starting points at a series of trial eccentricities, and set- 
ting the DM to the value measured earlier, we used TEMPO 
to compute xt surfaces (with 1^=22 degrees-of-freedom) as 
a function of P„i.b and x, by allowing the software to fit for 
/, the angle of periastron (oj) and the time since periastron 
(Tq). The values of Porb and x at the resulting xt minima are 
plotted in Fig.|6l and have typical errors in each parameter of 
10-20%. This process provided a 95% confidence lower limit 
to the eccentricity of e > 0.52. Orbits with e > 0.55 provided 
essentially perfect fits to the data with y^^jv ^ 1.0 and white- 
noise-like residuals with an RMS of ^30 /is. The spin-period 
behavior of four representative orbital solutions is plotted in 
Fig.lHas well as the measured spin-period behavior during the 
discovery observation. 

The family of orbital solutions that we have determined im- 
plies that M30B is part of a relativistic orbital system with a 
companion most likely in the mass range 0.25-1 M0. If the 
companion is a white dwarf (WD), the advance of periastron 
(w) due to general relativity will be several tenths of a de- 
gree per year and would be measurable (as will, therefore, the 
total system mass, M = m\ +m2) to a high degree of preci- 
sion within a year of monthly timing observations of a quality 
similar to the discovery observation. It is also possible, if 
such timing measurements can be made on a consistent basis, 
that additional post-Keplerian orbital parameters, such as the 
Einstein 7, will be measurable in M30B within several years. 
Such measurements would provide one of the few accurate 
mass determinations of a MSP which is important for test- 
ing models of NS spin-up via accretion. Most models predict 
the accumulation of >, 0.1M(7i of m aterial onto the NS dur- 
ing the spin-up process JPhinnev & Kulkarni 1994, and refer- 
ences therein), and recent observations have begun providing 
evidence for the case (e.g. Nice, Splaver, & Stairs 2003). 

In the current pulsar catalog, there are only four other GC 
binary pulsars with e > 0.1. The NS-NS system in M15, 
PSR B2127-11C (P = 30.5ms, P„rb = 0.335d, e = 0.681; 
iPrince etalJfT99ll) was almost certainly created in an ex- 
change interaction near the core of the cluster that replaced 
a lower-mass MS or WD companio n of a NS with another NS 
jPhinnev & Sigurdsson 1991; Sigurdsson & Phinnev 19931) . 
This interaction also placed the system as a whole in a highly 
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Fig. 6. — The family of possible Keplerian orbital solutions for M30B 
based on TEMPO fits to 26 TOAs determined during the discovery obser- 
vation. The 95% confidence lower-limit to the eccentricity is e > 0.52 (see 
j|5.H . ForPo,-;, (top) and asin(()/c (middle), typical uncertainties are 10—20% 
in the plotted value. In the bottom panel, the lower Hne shows the minimum 
companion mass (i = 90°) while the upper line represents the mass corre- 
sponding to the median inclination, (' = 60° , assuming random inclinations. 
In all cases we have assumed a pulsar mass of mi = 1.4Mq. 



eccentric orbit around the center-of-mass of Ml 5, which ex- 
plains its location near the outskirts of the cluster today. An 
exchange interaction also provides the most likely explana- 
tion for the eccentricity of PSR B 1802-07 in NGC 6539 
(P = 23.1 ms, Parb = 2.62d, e = 0.212; D 'Amico et al.tllT99l 
— a system possibly very similar to M30B — and a possible 
explanation for the eccentricity of PSR B1516+02B in M5 
(P = 7.95 ms, P orb = 6Md, e = 0.138; Anderson et al. 1997| 
as well. iRasio & Heggid ( 119951) have argued, however, that 
M5B's eccentricity may be the product of many distant in- 
teractions with other stars near the clusters core over several 
Gyr. The fourth cluster binary pulsar with a large eccentric- 
ity is the recently discovered PSR J 1750 -37 in NGC 644 1 
(P= 111.6ms, Porb= 17.3d, e = 0.71; Possenti et al.'"2001"), 
which may be a very similar system to M30B, albeit with a 
much slower spin period. Hypotheses of its origin await deter- 
mination of its position with respect to the cluster center, but 
with such a large eccentricity, an exchange encounter seems 
likely. 

We can estimate the time required for a pulsar with an 
initially near-circular orbit in the core of M30 to accumu- 
late an eccentricity of e > 0.5 by using Equation 5 from 
Rasio&Heggie (1991. With a one-dimensional veloc- 
ity dispersion for the cluster core of V£(0)=9.4±2.5kms"' 
tiGebhardt et al. 1995) . and the num ber density of s tars in the 
core of ~ 1.6 X 10^ pc"-' (tGuhatha kurta et alJll998h . we esti- 
mate ti,>o.5 — TPo^b^d Gyr, where Porb.d is the orbital period of 
M30B in days. Therefore, if M30B resides in the core of M30, 
it is possible that its eccentricity is due to many distant interac- 



Binary Millisecond Pulsars in M30 



9 



tions. It is interesting to note that for M30A, which probably 
does reside in the core of M30, the time required for multi- 
ple distant interactions to produce an eccentricity that would 
be detectable in our current measurements (e ^ 2 x 10""^) is 
~2-3Gyr. This means that M30A is either younger than 
~2-3 Gyr, or the mechanism described by Rasio & Hegg ie. 
(0.995) is not as efficient in M30 as we calculate. 

Should M30B reside outside of the cluster core, that may 
indicate that its high eccentricity was caused by an ex- 
change interaction similar to that which could have produced 
PSR B 1802-07 and almost certainly did produce the recently 
discovered eclips ing MSP in the outskirts of NGC 6397, 
PSR J1740-5340 (iD'Amico et al.l2001b) . If the exchange in- 
teraction replaced a low-mass WD with a MS star, optical ob- 
servations should allow its identification, once a timing posi- 
tion has been measured. Such a MS or even post-MS compan- 
ion could help to explain our lack of detections of the pulsar 
as well. If the companion star experiences significant though 
possibly erratic mass loss, the pulsed flux from M30B could 
be eclipsed or even quenched entirely for periods at a time, 
just as observed from PSR J1740-5340. If this is the case, 
higher frequency observations of the cluster may allow us to 
peer through the wind of the companion and see the pulsar 
consistently. 

If the companion to M30B is a MS or post-MS star, though, 
orbital circularization due to a hydrodynamical mechanism 
(Ijassoul 1995) may occur on timescales much shorter than the 
^j^Qg-iOyj. j.jjg system we would expect if it was cre- 

ated via an exchange interaction. Tassoul ( 1995) calibrated 
a theoretical relation for hydrodynamical orbit circulariza- 
tion by comparing the shortest orbital periods of MS bina- 
ries with eccentric orbits in several stellar clusters of different 
ages. He determined the time for circularization in years to be 

roughly f„>c ~ 1.5 x 10^^~'^/'*-*-P^^/^^, where the fitting param- 
eter ~ 8.3 for stellar systems of age >1 Gyr and P„rb,d is 
the orbital period in days. For the family of orbital solutions 
we have found for M30B, this corresponds to tdrc ~ 10^"^ yr, 
and implies either that the companion star is a WD, or that 
the interaction that exchanged a MS companion into the sys- 
tem happened within the last 10^ yrs, or significantly more re- 
cently if the orbital period of M30B is <6 d. We note, though, 
that the Tassoul (1995) calculations assume a binary system 
of two MS stars and almost certainly need modification for 
NS-MS systems. 



6. CONCLUSIONS 

We have used the Green Bank Telescope to discover two 
binary MSPs in the core-collapsed GC M30. One of these 
systems is a member of the rapidly growing class of eclips- 
ing MSPs and the other seems to be in a highly eccentric 
and relativistic orbit. Higher sensitivity observations of M30 
in the near future — for instance, using all of the available 
bandwidth provided by the GBT's 20-cm receiver — should 
allow us to monitor M30A for long-term variations in its or- 
bital parameters as has been seen in othe r eclipsing MSPs (e.g. 
lArzoumanian. Fruchter. & Tavloill9'9 4!) and enable us to con- 
sistently detect and time M30B. Additional observations at 
other radio frequencies will allow us to probe the eclipse re- 
gion of M30A to a greater degree, and will hopefully allow us 
to constrain the mechanism behind the eclipses themselves. 

It seems clear that as in the case of most other GCs (e.g. 
Camilo et al. 2000), we are only seeing the most luminous 
pulsars contained in M30. Given the extreme scintillation 
events we have witnessed from M30A and B, we consider it 
quite likely that the improved observations mentioned above 
will uncover several additional pulsars in M30 in the years 
ahead. 
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